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ABSTRACT 

' Aims. RR Lyrae stars have been observed to improve the insight into the processes at work in their atmospheres. 

' Methods. Simultaneous Stromgren-photometry allows to obtain a rapid sequence of measurements in which photometric indices are 

unaffected by non-optimum observing conditions. The indices y, b-y, and cj are used with an established calibration to derive T^if and 
to derive the gravity, log ggj from the Balmer jump, throughout the pulsation cycle. By employing the equations for stellar structure, 
additional parameters can be derived. 

Results. Stromgren photometry and its calibration in terms of T^g and log g can be used to determine the run of R and the atmosphere 
\l , pulsation velocity. We find that the Balmer-line strengths are correlated with T^ff and that the strength of the CaiiK line correlates 

. well with the radius of the star and thus the pulsation-dependent density of the atmosphere. The density in the stellar atmosphere 

fluctuates as indicated by the changes in the gravity loggBj, derived from ci, between 2.3 and 4.5 dex. Also the Stromgren metal 
index, mi, fluctuates. We find a disagreement between logg(r, L, M), the gravity calculated from T^g, L, and the mass M, and the 
gravity loggBj- This can be used to reassess the mass and the absolute magnitude of an individual star The curves derived for the 
^0 ' pulsational velocity Vp„i differ from curves obtained from spectra needed to apply the Baade-Wesselink method; we think these 

^' 



differences are due to phase dependent differences in the optical depth levels sampled in continuum photometry and in spectroscopy. 
We find an atmospheric oscillation in these fundamental mode RR Lyrae stars of periodicity P/l . 

Conclusions. Carefully conducted Stromgren-photometry allows to derive a large number of parameters for RR Lyrae stars. It provides 
Q . a means of deriving masses and absolute magnitudes. When comparing photometry results with spectroscopic analyses it appears that 

i-^ • optical depth effects affect all interpretations. 

52 ' Key words. Stars; atmospheres - Stars; variables, RR Lyrae - Stars; horizontal-branch - Techniques; photometric - Radiative transfer 

^ ' Stars; individual; RR Gem, SZ Gem, SY Ari, AS Cnc, BH Aur, TZ Aur, TW Lyn, CI And, AA CMi, AR Per, X CMi, BR Tau 

1^ 1. Introduction Using photometry of RR Lyr stars, one can derive the varia- 

' tion in the parameters of the stellar surface. The first investiga- 

^ ; When lower mass stars evolve and come to core He burning, tions of this kind were carried out by Oke & Bonsack (1960), 

^ their properties can become those of the "instability strip" in the Qke etal. (1962), Oke (1966), and Danziger & Oke (1967), 

. HRD. The atmospheric structure then reaches, as it were, a cer- who used spectral scanner data in comparison with model atmo- 

VO tain "undecidedness", causing rhythmic attempts to reach a sta- spheres to obtain T^g and log g to calculate the change in radius, 

ble state. The atmosphere cannot find its "thermal equilibrium" They noted that the changes in the atmosphere imply that one 

in the sense explored at a general level by Renzini et al. (1992) observes light from diff^erent atmospheric layers so that the radii 

in the gedankenexperiment the "gravothermal hysteresis cycle", jjjus derived might not be reliable. 

The rhythmic expansion and contraction of the atmosphere The information can also be obtained from Stromgren pho- 

of pulsating RR Lyrae stars is caused by the /^-effect. The optical tometry. After the establishment of the Stromgren-photometric 

5^1 depth K in the layer in which He is ionized has a rhythmic varia- system and its early calibration (see, e.g., Breger 1974) this pho- 

C3 ^ tion: a higher level of ionization causes lower opacity, leading to tometry was used by van Albada & de Boer (1975) to derive the 

a higher photon throughput, causing local cooling, which leads pai-ameters = 5040/reff, log g, and R for all phases of the pul- 

to recombination and thus to higher opacity, reduced radial en- sation. McNamara & Feltz (1977) and Siegel (1982) also used 

ergy transport so to increasing inner temperatures, coming full the Stromgren-system. Similai- studies were performed by Lub 

circle to increase ionization (see, e.g.. Cox 1980; Gautschy & (1977a, 1977b, 1979) with the Walraven photometric system. 

Saio 1995). This cyclic behaviour produces hysteresis effects in 1 f c nm . i. »u 1 f ■ 

. ,. ■' ^ , ^ „ . „ . Only for a few RR Lyrae stars has the cycle of variation 

the colour indices of the emergent stellar radiation as well as in , ■ • u » » • j » -i • »u c/ •• 

, ,,. , N r^, been followed in photometric detail in the Stromgren system, 

vanous other observables (such as spectral line strengths). These rr, , f f . • . . °. 

^ , ^. . ^ • , Thus only tor a tew of these stars is the variation in parame- 

all are based on, of course, hysteresis in surface gravity and ef- ^ , , .11 nu . . • 

fective tem r r = ter values over the cycle accurately known. Photometry is nor- 

ective temperature. mally performed sequentially in whichever wavelength bands se- 

lected. Sequentiality mandates good photometric conditions to 

* Original data (photometric indices) for each star are available in obtain accurate colour indices for the derivation of reliable stel- 

electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr lar parameters. Moreover, sequential measurements in the cho- 

(130.79.128.5) or via http;//cdsweb.u-strasbg.fr/cgi-bin/qc at?J/A-KA/, sen bands are asynchronous among the bands. Performing mea- 
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Fig. 1. Lightcurves of the RR Lyrae stars in y. Data from sub- 
sequent cycles have been combined (see Table[T]for observing 
dates). The curves have been shifted arbitrarily in y to fit them 
into one panel. For the stars lower in the figure, the lightcurve 
coverage is incomplete. 
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Fig. 2. Lightcurves of RR Gem and SY Ari observed with Busca 
in the Stromgren-system. Note the different brightness amplitude 
of the two stars. At the bottom, the instrumental brightness of 
the mean of two comparison stars (selected in the field of each 
variable) is shown. 



surements simultaneously in well selected photometric bands al- 
lows to obtain precise colour indices even in poor photometric 
conditions. For the research presented here, we used the Bonn 
University Simultaneous CAmera, Busca (Reif etal. 1999). To 
transform the colour indices to astrophysical parameters, the cal- 
ibration of the Stromgren system by Clem et al. (2004) has been 
used. Busca also allows a rapid succession of measurements 
leading to a better coverage of fine structure in light curves. 

The Baade-Wesselink method (Baade 1926, Wesselink 
1946) requires the measurement of the radial velocity variation 
to produce a full characterization of a pulsating star. Following 
van Albada & de Boer (1975) we will again show that, with ac- 
curate photometry, R can be derived (rather, the apparent angu- 
lar extent, which indicates the radius when the distance to the 
star is known) and the changes in R can be used to calculate the 
run of atmospheric velocities Vpui through the pulsational cycle. 
Stromgren photometry allows, in principle, to derive astrophys- 
ical parameters more accurately than possible from the widely 
used Johnson photometry. 

2. Observations and data reduction 

The observations were carried out in three sessions at the 
Calar Alto Observatory of the Deutsch-Spanisches Astronomic 
Zentrum (DSAZ) in southern Spain in three successive win- 
ter seasons, between 4 and 9 Jan. 2005, between 16 and 19 
Dec. 2005, and between 26 and 28 Nov. 2006. 

The times of maximum brightness within the observing ses- 
sions were calculated for each star from the period P and a refer- 
ence epoch as Vm^x - JD(ref.epoch)) x (nP) with appropriate n. 
This allowed optimal phase-coverage planning of the photome- 
try. In the various figures to follow, the data are plotted against 
phase <1) based on the observed maximum. For a few stars, the 



maximum could not be measured so the phase was determined 
from early epoch maxima. 

Our programme consisted of observing 18 relatively bright 
RR Lyrae (for accurate positions and additional data, see Maintz 
2005). Of 12 of these, (almost) complete lightcurves were ob- 
tained (coverage >70%) as given in Table[T] and Fig.[T] (and 
Fig.|2]i. However, for X CMi we missed the maximum. Of 6 stars 
(GMAnd, OVAnd, GIGem, TYCam, TTLyn, CZLac), only 
smaller portions (15 to 57%) of the light curve could be obtained; 
these stars are not discussed further. 

The stars selected are fundamental mode pulsators (see 
Fig EJ. They generally have lightcurve bumps. 

Busca (Reif etal. 1999) is operated at the 2.2m telescope 
of the Calar Alto Observatory. Busca splits the telescope beam 
above the Cassegrain focus via dichroic beam splitters in 4 wave- 
length channels between 3200 and 9000 A, with edges at 4300, 
5400 and 7300 A. Each channel is equipped with a (4kf CCD. 
In each channel, a filter wheel allows the placement of appropri- 
ate filters. 

For the measurements we used the Stromgren u,v,b,y- 
system, extended by the Cousins / bancfl. The y band is close 
to the cut-off of a BuscA-dichroic filter; however, the transfor- 
mation of y into the standard system could be performed with- 
out problems. Because of the proximity of the wavelengths of 
u and V, their filters are in the same Busca channel. Since the u 



' Since the Cousins / is about 10 times wider than the bands of the 
Stromgren-system, the / measurements were mostly overexposed and 
have not been used for this paper. For a few stars we also made some 
measurements in the H-Balmer filters H/3W (wide) and H/3N (narrow), 
which are in the same channel as the Stromgren-fo filter. Since the H/JN 
(narrow) filter is yet much narrower than the Stromgren-filters, exposure 
times would have to be excessively long and coverage of the light curves 
in these filters poor so these data have not been used either. 
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Table 1. The stars with basic information, the J.D. of the observations, and the phase ranges covered 



Star 


CI And 


SY Ari 


ARPer 


BRTau 


BH Aur 


TZ Aur 


RA 


01 55 08.29 


02 17 34.04 


04 17 17.19 


04 34 42.89 


05 12 04.26 


07 11 35.01 


DEC 


43 45 56.47 


21 42 59.26 


47 24 00.63 


21 46 21.72 


33 57 46.95 


40 46 37.13 


mean y ; E{B - V) [mag] 


11.86 ; 0.22 


12.21 ; 0.25 


9.40 ; 0.73 


11.96; 0.31 


11.43 ;0.19 


11.43 ; 0.19 


distance [pc] ; [M/H] 


1741 ; -0.83 


2100 ; -1.4 


612 ; -0.43 


1835 ; -0.7 


1400 ; +0.14 


1500 ; -0.79 


period [d] 


0.484728 


0.5666815 


0.425551 


0.3905928 


0.456089 


0.39167479 


exp. time [s] v, b,y; u,b,y 


40; 40 


70; 70 


20; 60 


70; 210 


70; 210 


30; 30 


start of observations 














at JD = 2453370. 


+5.2914 


+5.2947 








+5.6153 


phase covered 


0.82 to 1.26 
+7.4215 
1.21 to 1.43 


0.66 to 1.05 
+7.2612 

0.13 to 0.65 
+8.3353 

1.03 to 1.16 

6S to 91 








0.67 to 0.95 
+6.4747 

0.86 to 1.60 
+7.5535 

0.62 to 0.97 

1 g T 1 n/: 
-f-o.-7 lUO 

S5 tn 89 




+ 1.32 






j.^ ?QfiQ 


+1 9QSn 

^ 1 .Z.70U 




pilaSC LOVCICU 


n fifi tn n 






n S'^ tn 1 14 


n "^7 tn 1 
U.JZ LU L.jy 

+9 SRfiR 

^4 tn S4 

U.^T^ LU U.^T^ 










+7.2981 


+7 S01 ^ 

^ / . JU i J 






UiiaaC LUVCiCQ 






75 fn 1 fi^ 


n 7^ tn n 87 

\J. J J LU U. O / 






Star 


AACMi 


RRGem 


XCMi 


TWLyn 


SZ Gem 


AS Cnc 


RA 


07 17 19.17 


07 21 33.53 


07 21 44.62 


07 45 06.29 


07 53 43.45 


08 25 42.11 


DEC 


01 43 40.06 


30 52 59.45 


02 21 26.30 


43 06 41.56 


19 16 23.93 


25 43 08.80 


mean y ; E{B - V) [mag] 


11.13 ; 0.16 


11.03 ;0.22 


12.23 ; 0.35 


11.91 ; 0.16 


11.45 ;0.13 


12.58 ; 0.16 


distance [pc] ; [M/H] 


1220 ; -0.15 


1260 ; -0.29 


1766 ; -0.71 


1549 ; -0.66 


1601 ; -1.46 


2324 ; -1.89 


period [d] 


0.476327 


0.397292 


0.57138 


0.481862 


0.5011270 


0.61752 


exp. time [s] v,b,y\ u, b, y 


40; 40/20; 60 


20; 20 


70; 70 


40; 40 


20; 20 


40; 40 


start of observations 














at JD = 2453370. 




+5.5257 


+8.6759 


+6.5937 


+5.5385 


+7.5695 


phase covered 




0.98 to 1.49 


1.08 to 1.17 


1.16 to 1.51 


0.39 to 0.77 


0.73 to 1.05 






+6.4776 


+9.3641 


+7.2652 


+6.4800 


+8.4688 






0.38 to 1.10 


0.29 to 0.93 


0.55 to 0.76 


1.26 to 1.48 


1.18 to 1.63 






+8.5121 




+8.3421 


+7.6349 


+9.4268 






0.50 to 0.88 




0.79 to 1.25 
+9.2625 
0.70 to 0.76 


0.57 to 0.84 


0.73 to 1.28 


at JD = 2453720. 


+2.4285 




+ 1.5757 




+ 1.5348 


+ 1.4625 


phase covered 


0.77 to 1.09 




1.21 to 1.32 




0.82 to 1.28 


0.59 to 0.69 


at JD = 2454060. 


+6.6066 




+6.6110 








phase covered 


1.02 to 1.32 




1.07 to 1.28 









[M/H] (in logarithmic units) estimated from our nii (see Sect.O. 




0.5 - 



0.2 0.4 0.6 0.8 
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Fig. 3. Plot of the amplitude in y versus period of our stars show- 
ing that they are fundamental mode RRa-type pulsators. The 
band between the dotted lines marks the zone of RRa and RRb 
pulsators (see Ledoux & Walraven 1958). 



and V bands could not be measured simultaneously, the photom- 
etry alternated between the simultaneous exposures in y, b, v, (/) 
and y, b, u, (/). For two stars, light curves in the four Stromgren- 
bands are shown in Fig.|2] More data are shown in Maintz 
(2008a). 



To calibrate the data, the instrumental brightness of each 
variable was obtained by comparison with two non-variable stars 
in the field. Where possible, one star with a brightness near the 
Jmax of the RR Lyrae star was used and another of brightness 
close to jmin of the variable. (In the data a known suspected vari- 
able was quantified; see Maintz 2008b.) The airmass correction 
was obtained from these comparison stars. 

To secure the absolute calibration, we observed appropriate 
stars from the list of Perry etal. (1987) at moments in the light 
cycle of the variable where the changes in brightness were small 
or monotonous. The instrumental magnitudes of the (airmass- 
corrected) photometry were related to the true magnitudes of the 
calibration stars. This permitted the calibration of the RR Lyrae 
star photometry. First, the RR Lyrae y measurements were cal- 
ibrated. Since the instrumental (b - 3')inst from Busca is more 
accurate than bcM-ycM (the same is true for the other indices), the 
calibrated b - y and v - b were obtained from the instrumental 
index and that of the calibration stars. The sequentially obtained 
colour indices u - y and v - y were then used to interpolate in 
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Fig. 4. Lightcurves in the calibrated Stromgren-indices of the 
stars RR Gem, SY Ari, and BR Tau. The phase coverage of the 
photometry can be reconstructed from the information given in 
Table[T](for BR Tau the coverage is only 70%) and can be recog- 
nised in Fig.[T] 



time to obtain u - v, which was calibrated. Using these indices, 
ci = (m - v) - (v - b) and nii - {v-b)-{b- y) were calculated. 

The RR Lyrae star CCD photometry was performed with 
equal exposure times throughout the cycle, of 20 s for the bright- 
est and 70 s for the faintest stars. Exposure times are given in 
Table[T] For three stars the exposure times in the u, b, y measure- 
ments were 3x those of the v,b,y measurements (see Table[T]). 
The (simultaneous) reading-out of the CCDs took about 2 min. 
Thus the CCD exposures followed each other quickly, within 
2.5 min for the brighter stars and 3.5 min for the fainter ones. 
The short exposure times did not lead to large time smearing, 
not even when the star changed quickly in brightness. In most 
cases, two stars were observed alternatingly in rapid succession 
to achieve good light curve coverage for as many stars as pos- 
sible. This alternation explains why most light curves presented 
have data intervals longer than 3.5 min. 

Basic data {d, E{B - V), [M/H]) were taken from Beers et al. 
(2000), Fernley etal. (1998), and Layden (1994), in part taking 
averages. Some of the stars are known to be behind gas with 
extinction and we corrected for that using expressions given by 
Clem et al. (2004). For four stars we calculated d from My based 
on [Fe/H] using the formula of Fernley et al. (1998). The adopted 
distances are given in Table[T] 

3. Colour-index loops 

For three stars the run through the cycle of the derived absolutely 
calibrated colour indices is shown in Fig.|4] The b - y versus 
c\ diagram for two typical stars from our sample, RR Gem and 
TZ Aur (Fig.|5] left panel), shows clear colour index loops. The 
index mi varies in line with c\ (FigHJl. Figure|5] (right panel) 
shows the very moderate m\ loops in a colour-colour plot. 

The {ci,b-y) curves show the interplay of the surface param- 
eters \ogg and Teff during the cycle. The loop is essentially due 
to ci , which represents (above T^s ^ 6000 K) the strength of the 



Balmer jump. At low Teff, the Balmer level is hardly populated. 
When the temperature increases (bluer b - y), T^s and (and 
thus log g) of the atmosphere gas increase and (collisional) ex- 
citation of hydrogen becomes noticeable. When the temperature 
decreases the population of the Balmer level decreases as well. 
In RR Lyrae star atmospheres, the ionisation of hydrogen hardly 
plays a role (except for an increase in Ug) because the highest 
temperatures in the pulsational cycle are Teff < 9000 K. 

It is evident from the data that the highest temperatures ap- 
pear when the star is close to maximum light. During most of 
the period, the measured colour indices indicate low tempera- 
tures (large b-y, small c\ ), as visible in the clump of data points 
in the lower-right corner of the left-side panels of Fig.|5] The 
brightening leg and the dimming leg of the colour curves do not 
overlap. During brightening, b-y changes faster than ci, which 
is indicative of a faster rise in Teff than in the excitation to the 
Balmer level. In the dimming branch the converse is found (see 
Fig.|5] left). Thus the colour curve exhibits hysteresis. 

As for the loops in m\ (Fig.|5]right), these are quite moderate 
and do not show much in the sense of hysteresis. The change 
in m\ over the cycle is on average 0.1 mag (;«i is largest at 
maximum light). For our stars, nii is larger for stars known to 
have larger [M/H] (see also the calibration for red giant stars by 
Hilker 2000), while for stars with [M/H] < -1 dex mi is only 
marginally metal dependent. At the T^s of RR Lyrae stars, cycle 
variations in m\ related to the [M/H] are small. Nevertheless, us- 
ing these trends, we estimated [M/H] for the two stars for which 
these values were unavailable from other sources (see TablelTJ. 
The variation in mi over the cycle must be caused by changes in 
temperature and electron density, which both influence the level 
of ionisation and excitation of ions. 



4. Behaviour of atmospheric parameters 

4.1. Deriving parameters in tlie stationary case 

The parameters Teff and \ogg were determined from our pho- 
tometry using the conversion grids given by Clem etal. (2004). 
These grids are based on modelled spectral energy distributions 
and are provided in steps of 0.5 dex of metal content [M/H]. 
Clem et al. made extensive calibrations of the Stromgren system 
against the parameters of non-variable stars. They show that in- 
dividual Population II star colours can be matched with models 
to within ^200 K and ^0.1 dex, showing what is the ultimate 
uncertainty of a match. For Population I stars, the effect of small 
variations in metal content (0.3 dex) on the colours, and thus on 
the derived Tgff and log g, is small. The dependence of colours 
on metal content is significant only for Teff < 6000 K. For each 
star we used the grid for which [M/H] was closest to the metal 
content of each star. 

The parameters Teff and \ogg can be used to derive ad- 
ditional parameters describing the behaviour of the stellar at- 
mosphere. The conditions in a regularly pulsating stars are in 
"quasi-hydrostatic equilibrium" (Cox 1980; Ch. 8.3), meaning 
that the use of a calibration based on stable stars is justified. 

However, when a RR Lyrae star brightens from Vmin to Vmax 
within approximately one hour, the atmosphere is probably not 
in quasi-hydrostatic equilibrium. Doubling of Ca and H lines has 
been seen (e.g., Struve 1947, Sanford 1949) and H lines have 
been seen in emission (e.g., Preston & Paczynski 1964) in that 
part of the RR Lyrae cycle, which according to Abt (1959) origi- 
nate from dynamic effects. Since no models exist that take these 
effects for photometry into account there is no other option than 
to use static models. Furthermore, the light we detect comes 
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Fig. 5. Colour loops in b - y,ci (at left) and in b - y,mi (at right) of several stars. The data points of the brightening phase are 
marked with •. 

Left: Clear hysteresis-like behaviour in b - y, ci is shown for RR Gem and TZ Aur. Arrows mark the direction of time. At phases 
between 0.3 and 0.85 (see Fig.|4|i, the colour index ci seems to just scatter The b-y,ci combination represents mostly Teft but also 
Balmer excitation. 

Right: Colour loops in b - y, nii (shown for RR Gem, AS Cnc, SY Aur, TW Lyn, and SZ Gem) marginally exhibit hysteresis-like 
behaviour The index m\ generally represents metal content [M/H], its change in RR Lyrae stars is mostly due to changing T^g. 



from the photosphere, i.e., the layer in which, at the particular 
wavelength observed, the optical depth is t ^ 0.7. This means 
that in the course of the pulsational cycle it is not necessarily the 
same gas from which the light detected emanates. We return to 
these aspects later (Sect. 18. 2]) . 

The luminosity of a star is related to the surface parameters 
R and Teff through the familiar equation L - AnR^ ■ crT^^, which 
can be rewritten in the relative logarithmic form 



(1) 



The surface gravity of a star is defined as ^ = G M/R^. This 
expression can be transformed into a logarithmic one relative to 
solar values given by 



(2) 



Using Eqs.[T]and|2l one can then eliminate the radius so that 
the combined equation has the stellar mass, temperature, gravity, 
and luminosity as variables, 



log (^) + log g = 4 log Teff + log - 10.68 



(3) 



where -10.68 - (logTeff - logg)Q in cgs units. This equation 
can be used to calculate log g. 

4.2. Deriving parameters in the puisating case 

In a pulsating atmosphere the parameters Tetf , log g, and R vary 
continuously (note that no time subscripts are used in this paper). 



• Teff, L. The time-dependent values L and Teff can be derived 
from the photometry if one knows the stars distance. Teff is sim- 
ply inferred from a calibrated b-y (using the grid of Clem et al. 
2004). L has to be determined from y, the integral over the spec- 
tral energy distribution and the distance d of the star. 

As for L, an RR Lyrae star has the convenient property that 
the maximum of its spectral energy distribution (in B^) lies close 
to the middle of the visual wavelength band, near Johnson V or 
Stromgren y. This is true during the entire cycle: Tetj actually 
varies only between roughly 5000 and 9000 K. This means that 
the bolometric correction is neither large nor varies a lot. We 
have performed the bolometric correction using the classic val- 
ues of Schmidt-Kaler (1982). We fitted a quadratic equation in 
the temperature range relevant for RR Lyrae stars. 

Distances were adopted as described in Sect.|2l If the dis- 
tance c/ of a star were off by 20%, this would result in an error of 
0.08 in logL. 

Results for L and T^ff are shown in Fig.|6] 

• R. The radius can be calculated from Eq.[T] For examples of 
radius change see Fig.|7] 

• log g, log gB], log ggg. The surface gravity g can be calculated 
from Eq.[3] if one knows the luminosity of the star (discussed 
above) and its mass (see below). This gravity is like an overall 
equilibrium gravity and is henceforth called g(T, L, M). 

There is a further aspect affecting the gravity. In a pulsating 
atmosphere, the actual surface gravity (g{T,L,M) at the level 
T ^ 0.7) is modified due to acceleration of the atmosphere during 
the pulsation cycle with respect to the "normal" gravity. Thus, g 
is modified by an acceleration term, (fiR/dt^. The total gravity is 
called the "effective" gravity 



?,s^g(T,L,M) + ^ 



(4) 
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Table 2. Cycle averages of several RR Lyr star parameters 



Star 


(L) 




log{g(T,L,M)) 


l0g<gBJ> 


RRGem 


43.5 


6647 


2.84 


3.33 


TWLyn 


31.6 


6196 


2.85 


3.61 


AS Cnc 


49.3 


6592 


2.77 


3.46 


SY Ari 


50.0 


6059 


2.61 


2.93 


SZ Gem 


52.9 


6739 


2.78 


3.38 


BH Am- 


46.4 


5997 


2.63 


3.04 


TZ Aur 


50.3 


5934 


2.58 


2.86 


AR Per" 


56.3 


6030 


2.55 


3.66 


CI And* 


47.7 


6270 


2.69 


3.19 


BR Tau'^ 


44. 


5980 


2.65 


2.92 


AA Cmi' 


45. 


6340 


2.73 


3.28 



Notes. L in Lg; r in K; g in cm s"^. L derived from y. Ay, d, and B.C. 

Light curve sparsely measured (see Fig.[T](. Light curve with a 
large gap in the descending branch; interpolated. Relatively large ob- 
serving gaps (see Fig.[T](; interpolation uncertain. 



indicating the vertical gravitational force in the atmosphere. 

The time-dependent Stromgren photometry measures the ac- 
tual condition of a stellar atmosphere. The gravity derived from 
b - y,ci is a function of gas density, visible in the excitation of 
the Balmer level and in ci representing the Balmer jump. We will 
call this gravity the "Balmer-jump" gravity, loggBj- This grav- 
ity is independent of log geff and instead represents the pressure 
of the gas sampled, i.e., f gas - f{T,p). 

We will return to the values derived for log g{T,L,M) and 
log^Bj in Sect.|9] 

• My, M. RR Lyrae star distances have thus far been derived 
using a reference value for My- This value of My does not apply 
to all RR Lyrae stars since horizontal-branch (HB) stars evolve 
over more than AlogL - 0.1 on the HB. If, e.g.. My of a star 
were to differ by 0.25 mag from the reference value adopted, L 
would be off by 0.10. 

The mass of RR Lyrae stars is in the range from 0.6 to 
0.8 Mq. We adopted a mass of M = 0.7 Mq, thereby considering 
that our RR Lyrae stars can deviate 15% from that value. This 
deviation propagates to a maximum deviation of 0.06 dex in a 
calculated log g{T, L, M). 

A few observational determinations of the absolute magni- 
tude My and mass M of horizontal-branch and RR Lyrae stars 
exist, such as those of de Boer etal. (1995), Moehler etal. 
(1995), de Boer etal. (1997), Moehler etal. (1997), Tsujimoto 
etal. (1998), and Gratton (1998). 

We return to the effect of incorrect values of these reference 
parameters in Sect. 14. 51 

The run of T^s, log g^], and R derived as described above are 
shown for three of our stars in Fig.|2] The shape of these curves is 
very similar to those presented by van Albada & de Boer (1975). 
But the data of Fig.|7]is less noisy because the photometry with 
CCDs is faster and more precise than possible with the slower 
photomultipliers of that time, and because of the simultaneity of 
the BuscA measurements. 

4.3. Phase resampling in steps of 0.02 

To facilitate the calculation of averages over the phase and addi- 
tional analyses, we resampled the curves of L, Teff , R, and log ^bj 
as derived from the observational data using steps of AO = 0.02 
(starting with / - near maximum light). These resampled val- 
ues (given henceforth with subscript 0.02) can be easily inte- 



grated (fixed time step). The result of these resamplings can be 
seen as curves in Figs.|6]and|2](and in further figures). 

Time averages over the cycle of our stars have been calcu- 
lated based on the 0.02 phase stepped lightcurves. For a few 
stars, the lightcurve coverage was not complete and we made 
reasonable interpolations in the gaps to obtain the time averages. 

4.4. Time-averaged quantities 

The structure of a "non-pulsating RR Lyr star" is given by the 
time-averaged parameters (Tetf), (L) and (g). The time-averaged 
gravity, {g), has in a purely periodic system zero contribution 
from d^R/dt^. 

Oke and collaborators and van Albada & de Boer felt that 
since g varies strongly, the "average" value is best derived from 
the quiet part of the cycle, i.e., from 0.3 < O < 0.85. However, 
Fig. |9] makes clear that this has shortcomings. In that phase in- 
terval, log ^(r,L,M) is (for most stars) lower than average by 
between 0.1 and 0.2 dex, because the atmosphere is then slowly 
but continuously contracting (but see Sect.|9]). 

The time-averaged quantities of essential parameters were 
calculated from the resampled curves in steps of AO = 0.02 (see 
above), the integral was derived as the sum of 50 values. Thus 

{T,ff} = J Teff dt = (l^f J ( (reff)o.02) ,) / 50 (5) 

and 

<L) = (4°i (Lo.o2),)/50 , (6) 

individual values of L0.02 being calculated as described in 
Sect. l4.2l The averages of the two gravity versions are 

{g(T, L, M)) = (l^^'i ( g{T, L, M)o.o2) /) / 50 (7) 
and 

<gBJ> = (2-'\((^Bj)0.02),)/ 50 . (8) 

Time-averaged values for 1 1 stars are given in Table|2] and 
are for four stars included in Fig.|6] It can immediately be seen 
that the values of {g(T, L, M)) and (^bj) are quite dissimilar We 
discuss this further in Sect.|9l 

4.5. Error budget 

In our analyses, the following uncertainties had to be taken into 
account. 

Photometric errors are, given the instrumental characteristics 
of BuscA, small. We estimated the uncertainties in j to be 0.03 
mag. The bolometric corrections are small and introduce errors 
of up to 2%. This leads to a total typical photometric uncertainty 
in L of 3%. In the colour indices b - y, m\, and c\, the uncer- 
tainties are 0.02 mag (smaller than the error in y because of the 
simultaneity of the measurements). 

We used distances from the literature. The effect of their un- 
certainties on L is given below. 

To obtain Tetj and log ^bj, the indices b-y and c\ were used 
with the grid of Clem et al. (2004). The grid is not perfect and 
reading T^s and log ^bj from that grid by eye, we assessed that 
(based on readings performed three times for the entire cycle 
for one star) the uncertainty in Tgff is about 1 %, and in log ^bj 
about 0.05 dex. Since the effects of metallicity, [M/H], are only 
significant for T^s < 6000 K, which occurs for only 3 of our less 
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Fig. 6. For four RR Lyrae (TW Lyn, RR Gem, AS Cnc, and SY 
Ari), the run of Tetf and L is shown in an HRD. The originally 
derived values are plotted as 3-pointed stars, the resampled curve 
(in steps of 0.02 in phase, Sect. l4.3l l is shown as a full line. 
The mean values over the cycle (from Table|2]) are given as •. 
The background data (ZAMS, metal-poor ZAHB, the evolution 
tracks of a 0.55 and a 0.6 Mq HB star) are taken from de Boer 
& Seggewiss (2008; Fig 10.9). Two lines of constant R sue in- 
dicated. The almost vertical dashed lines mark the approximate 
location of the pulsational instability strip from Gautschi & Saio 
(1995). The stellar surface parameters of the RR Lyrae show hys- 
teresis. 

well-observed stars near minimum light, we are confident that 
our error estimates are only at those temperatures affected by a 
possibly small mismatch in [M/H]. 

Given the measurement uncertainties in L (3%) and Teff 
(1%), it follows that the uncertainty in the calculated R values 
is 4%. If, in addition, a distance were incorrect by H-10% that 
would give an additional offset for a given star in /? of -1-5%. 

As mentioned in Sect. l4.2l we assumed a fixed mass for all 
stars of 0.7 Mq, and indicated a margin of error of 15%. This 
error would lead for a given star to an offset of 0.06 dex in a 
calculated log g{T, L, M). 

The values of log g(r,L,M) calculated from Eq.|3]have a 
measurement error of about 8% or 0.03 dex. A distance error 
of H-10% leads to an offset in log g{T, L, M) of -0.04, an error in 
the mass M of H-10% to an offset in log g{T, L, M) of +0.04. 

Finally we note that the distances of RR Lyrae have es- 
sentially all been derived from a reference My (distance from 
distance modulus), in some cases adjusted for metal content, 
[Fe/H] . For a review of the latest calibrations we refer to Sandage 
& Tamman (2006). However, in spite of all calibration efforts. 
My has a spread in My of between 0.2 and 0.4 mag for a given 
[Fe/H] or an uncertainty of even up to 1.0 mag if [Fe/H] is un- 
known. Moreover, given the evolution of RR Lyrae stars on the 
HB, a single and universal value of My does not exist anyway; 
the evolution of an HB star may lead to an increase in log L of 
0.5 dex when reaching the terminal-age HB (see, e.g., de Boer 
& Seggewiss 2008; their figure 10.9). Therefore, if a considered 
star were brighter than the chosen reference value, its derived 
distance modulus and thus distance would be too small. If the 
true My of the star is unbeknown to be 0.25 mag brighter than 
the reference value, the distance (as given in Table[T]i will turn 
out to be too small by 12%. This affects the derived L, then be- 
ing too small by 0.10 dex. 
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Fig. 7. Run of the parameters Teff, log^Bj, and R (derived as de- 
scribed in Sect. 14. 21 ) during the pulsational cycle of RR Gem, 
AS Cnc, and TW Lyn using the calibration grid by Clem et al. 
(2004). The thin dashed lines are the values resampled at AO = 
0.02 (see Sect.l43]). 

The mean values of Table|2]were derived from a large num- 
ber of individual measurements, so the observational scatter is 
basically cancelled. What remains is the error related to the cal- 
ibration grid and also the error in the assumed distance and/or 
My (affecting log §bj and L as just indicated, as well as R). 

5. Hysteresis of the parameters Jeff, R, and loggBj 

The values of the atmospheric parameters for all points in the 
cycle were derived from the photometry as described above. In 
Fig.|6]we show the run of T^^ and L of four of our stars in an 
HRD. These parameters exhibit loops through the diagram ex- 
tending over and beyond the instability strip. The loopings are 
reminiscent of hysteresis. In the older literature, the lagging be- 
hind in the change of some parameter compared to another pa- 
rameter is called "phase-lag" (see, e.g.. Cox 1980). 

The behaviour of Tefr and loggBj against phase (Fig.|7]) 
shows remarkable features. For all stars, loggBj begins to 
increase about AO - 0.05 before Teff starts to increase. 
Furthermore, log §bj reaches its highest value when Teif has in- 
creased about halfway to its maximum value. And when T^s 
reaches its maximum, log^Bj is already decreasing (at phase 
^ 0). The resulting value of R (that of Fig.|7] is derived from 
Eq.[Tli indicates that R is smallest when log gBj is largest. 

The run of R shows that RR Lyrae are relatively extended 
stars for about half of the cycle (0.25 < O < 0.75), while large 
changes in radius take place between 0.75 < O < 1.25, with 
the most rapid change occurring between 0.9 < O < 1.15. At 
O ^ 0.9, a drastic reduction in the derived R takes place (the 
atmosphere apparently starts to collapse), reaching the smallest 
R near O ^ 0.95. At smallest R the atmosphere must have a high 
density. Up to the highest density, T^-s is at the same level as in 
the quiet phase leading up to the collapse. Soon after O ^ 0.95, 
R increases again. 
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In general, when y is brightest T^g is at its highest and R 
is very small; but immediately after minimum light, at the high 5 
peak in log gsu R is at its smallest. 



6. Spectral line strengths 

For three of our stars (TZ Aur, RR Gem, RS Boo), spectra were 
obtained (in January 2006, independently of the Stromgren- 
photometry) with the focal reducer spectrograph at the 1 m tele- 
scope of the Observatory "Hoher List" of the AlfA. Radial ve- 
locities could not be obtained with this instrument. The spectra 
have a resolution of ~3 A. Because of telescope size and focal- 
reducer set-up, exposure times had to be rather long, of the or- 
der of 3% of the period, leading to phase smearing. We deter- 
mined the equivalent withs of the lines of Ha, Hj8, Hy, Na i D, 
and Ca 11 K. 

For RR Gem and TZ Aur, we show (Fig. [8]) how the strengths 
of the Hy and CaK lines vary during the cycle. Results for 
RS Boo are not shown because the lightcurve coverage was poor 
Hy. The strength of Balmer lines is governed by T^g as ex- 
pected since the Balmer level can be populated only during the 
higher T^ff part of the RR Lyrae cycle. The temporal behaviour 
of W(Hy) and T^ff (Fig.[8]l is indeed very similar 
CaiiK. The strength of CaK is predominantly set by the gas 
density, i.e., the ionization balance forces Can to be reduced at 
higher density (when R is smallest). The curves of W(CaK) and 
R are very similar in terms of their temporal behaviour (clearly 
for RR Gem). 

The Na i D absorption is very weak and provides Uttle infor- 
mation, except that Na iD is strong when Tgff is at its lowest; the 
increased strength is caused by the shift of the ionization balance 
towards the neutral state at lower Tsg. 



7. The variation in the pulsation velocity Vpui 

7.1. The coarse variation in Vpui 

Figure]?] shows for three stars the change in radius, R, as derived 
from Eq.[T] One can thus calculate the velocity of the stellar at- 
mosphere (the "pulsation velocity", Vpui). We used 



Vpui = 



AR 
It 



R(b — Ra 



(9) 



where R is at intervals of fixed value AO - 0.02 (see Sect. l4.4b . 
To avoid being affected by peaks in the noise we applied a run- 
ning triangular (1,2,1) smoothing to the curves of R before cal- 
culating Vpui. In Fig.|9] Vpui (as well as the smoothed R) is given 
for eight stars with the best data. Note that Vpui plotted is that 
seen from the centre of the star (and is not a heliocentric value). 

The calculated variation in Vpui over the cycle (see Fig.|7]l 
has a general pattern. There is a pronounced minimum close to 
<5 ^ 0.95 where Vpui ^ -150 to -400 kms"', followed by a 
steep rise to a level of Vpui - -1-50 to 100 km s"' . In the course of 
the cycle, Vpui then slowly decreases and slowly becomes nega- 
tive reaching a level of Vpui - -50 kms"'. A sudden additional 
decrease sets then in to reach the most negative value of Vpui, the 
value this description started with. It shows that the atmosphere 
really collapses between O = 0.9 and 1 .0. Thereafter, the expan- 
sion is drastic but beyond O ^ 0.05 it is far more gradual (at 
Vpui - 50 kms"') as it decelerates. 

The behaviour of Vpui is similar to that of a model atmo- 
sphere as presented in the gedankenexperiment by Renzini et al. 
(1992). That gedankenexperiment had a star with a core and an 
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Fig. 8. For RR Gem and TZ Aur, the run of the equivalent widths 
W of the spectral lines Ca 11 K and Hy are shown and compared 
with the parameters R and Tgff (as in Fig.|7]i. The vertical scale 
is that of R [R©] of RR Gem. The values of the other parame- 
ters (W in A and T^ff in K) have been renormalized to fit this 
vertical scale (see the labels to the curves). For RRGem, Hy is 
proportional to Teff and is strongest when Teff is largest (stronger 
excitation); CaK is weakest when the radius R is smallest (at 
constant T and higher density, n^. shifts the ionization balance). 
The spectral observations of RR Gem come from two nights in 
the same time frame as but independent of the Stromgren pho- 
tometry (explaining the overlap of data in phase). In TZ Aur the 
relation of R with Ca K is less clear. The spectra of TZ Aur did 
not cover the entire light curve (gap near O = 0.15). 



envelope, in which the luminosity offered to the base of the en- 
velope. Lb, was manipulated. A gentle decrease in this Lb from 
a high Lb state could be accommodated by the envelope until 
some critical value of Lb at which the envelope collapsed. Then, 
in the gedankenexperiment. Lb was raised again and the now 
compact atmosphere slowly expanded until some other critical 
value of Lb was reached, now leading to a run-away expansion. 

In the case of RR Lyrae, the run of L is controlled by the 
opacity in the He^-He^^ ionization boundary layer With in- 
creasing transmission of energy through the envelope, L in- 
creases and the star expands until at some limit the envelope has 
become so tenuous that (since L can no longer increase) it must 
collapse. 

Having determined the run of Vpui one can derive These 
acceleration values change considerably during the pulsational 
cycle. However, because of the uncertainty in all Vpui values due 
to opacity effects to be discussed in Sect.[8j we refrain from com- 
menting on these acceleration values. 

7.2. The detailed variation in Vpui and a correlation with P 

The small-scale variations in the curves of Vpui (Fig-IS might 
be regarded as due to photometric noise. However, on closer in- 
spection there is significant structure. 

We selected the well-behaved Vpui-curves to determine an 
average Vpui-curve for an RRLyr star For that, we chose the 
stars RR Gem, TW Lyn, AS Cnc, and S Y Ari (the top row panels 
in Fig.|9|l. Before the curves can be added for averaging, small 
phase shifts were applied because the chosen y^ax was not per- 
fectly aligned with the epoch relation and the resampling at 0.02 
intervals in phase may also have introduced small phase shifts. 
Figure|9] shows the results for these stars. It is evident that the 
wiggles in the Vpui-curves are very similar for these four stars. 
For the four less well-observed stars not included in the averag- 
ing, a comparison with the average Vpui-curve shows that these 
stars have also the same behaviour in Vpui. 
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Fig. 9. The change in Vpui is shown together with the change in R (note that the scale of R differs from panel to panel). Top row: the 
"well behaved" and well-observed stars. Bottom row: the stars with less well behaved Vpui-curves (partly due to poor light-curve 
coverage and to interpolation). Top curve and data in each panel: The values of R were calculated from photometric data using 
Eqs.[T]and|2](as in Fig.|7]i. The resampled data (AO - 0.02) were used to derive the variation in R over these intervals. The curves 
shown give the resampled values after an additional running triangular (1,2,1) smoothing. Note that, for several stars, the period 
was not fully covered by our observations. Bottom curves in each panel: The run of Vpui as derived from the smoothed R-cnrve is 
shown. To guide the eye, a line at ypui=0 kms ' is added. The average of the Vpui curves of the four well observed stars (the stars 
shown in the top row) has also here been added as a dashed line (for each star a small appropriate phase shift was added). Clearly, 
the wiggles in Vpui are not noise but are systematic Vpui-variations. This is confirmed by the stars of the bottom row, whose Vpui 
curves were not included in the calculated average. The average Vpui curve shows a rhythm with a period of about P/1. 



For all stars, the wiggles in the average Vpui indicate that after 
the collapse the envelope oscillates (continues to oscillate) with 
a rhythm of about 1/7 of the period P. This rhythm at P/1 is 
the same for all our stars, independent of their actual period. It 
is as if the envelope has a residual oscillation triggered by the 
collapse. The collapse itself occurs from Vpui- through the 
highest collapse velocity back to Vpui- km s"', which is a time 
equal to about P/1. 

One might suspect that P/1 is due to our resampling at 0.02 
phase intervals, the period being (almost) exactly 7 of our data 
points. However, we note that Paparo et al. (2009) reported a spe- 
cific role of the 7th harmonic in one well-studied CoRoT RR Lyr 
star. The CoRoT data do not have the 0.02 phase resampling 
rhythm. Barcza (2003) found in the RRLyr star SUDra from 
photometry an "undulation"' of P/5. It is unclear what the cause 
of these periodicities is. 



8. Comparison with literature data and optical 
depth effects 

We derived the time-dependent radius of each star from our pho- 
tometry and from its change the velocity of the atmosphere. The 
analysis above was performed without considering analyses of 
other RR Lyrae data presented in the literature. 

A well-known method for deriving the radius of RR Lyrae 
stars is the Baade-Wesselink method, where the radial velocities 
observed over the cycle are integrated to infer the variation in ra- 
dius (see, e.g., Liu & Janes 1990). This method has been applied 
to several RR Lyrae stars. 



8.1. Comparison with otiier results 

For two of our stars (RR Gem and AR Per), cycle velocities were 
presented by Liu & Janes (1989). Their velocities were obtained 
by cross-correlating their RR Lyrae star spectra with spectra of 
two standard stars, one of spectral type F6 IV and one of A2 IV. 
The values they derived from the two standard stars are very 
similar so, that their final velocity curves are based on just the 
standard of type A2 IV. 

Our velocity curve for RR Gem has a shape quite different 
from the one of Liu & Janes (1989). Note that in their figure 5 
they give the observed radial velocity, whereas in our Fig. |9] the 
velocity as seen from the stellar centre is given. Their velocity 
curve shows a large speed away from the stellar centre at <1) ^ 
followed by a gradual slowing down to approach the minimum 
in velocity at O ^ 0.7. The velocity curves of other RR Lyrae 
stars of Liu & Janes show the same behaviour The amplitude in 
their velocity curves is ^60 km s"^ Adopting the usual geometry 
factor p = 1.32 (see Liu & Janes 1990; but see Fernley 1994) 
to convert from observed radial velocity to pulsational velocity, 
their velocities translate into an amplitude in Vpui of -80 km s"' . 

On the other hand, our stars have amongst themselves similar 
behaviour Their outer layers experience during the cycle only 
a modest change in velocity except for a rapid fall toward the 
stellar centre as of O ^ 0.85 to a fastest centre approach at O ^ 
0.95, followed by a rapid decline in velocity to reach a modest 
expansion velocity near O ^ 0.05. In our data, the full amplitude 
of the velocity Vpui is -300 km s"' . 

The difference between our velocity curves and those of Liu 
& Janes is mostly caused by the drastic reduction of R we find 
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Fig. 10. Run of the two aspects of log g (resampled data, see Sect. l4.3l l for the well observed stars. Top curves: The slowly varying 
curve is log g{T, L, M) derived from Eqs.[T] and |2] (thus using T and L derived from y,b - y), and the full horizontal line gives the 
average of {giT, L, M)) over the cycle (in logarithmic form; see Table|2]l. The peaked curve shows the gravity, loggBj, derived from 
the b - y,ci photometry and the grid calibrated to obtain Teff and log^Bj- Bottom curves: the ratio g{T, L, M)/gBs is shown on a 
logarithmic scale. Clearly, when A = log(g{T, L, M)/gBs) < the atmosphere gas must be quite compressed. The curve of Vpui has 
added (in units as in Fig.|9j. 



at O ^ 0.9. Eliminating this behaviour from our R curves (elim- 
inating the R data for 0.85 < O < 0.05) results in a velocity 
amplitude of Vpui of up to ^90 kms"', quite similar to the one 
of Liu & Janes mentioned above. In short, the radius change 
of an RR Lyrae star based on the Baade-Wesselink radial veloc- 
ity curves does not chime with the pronounced change in R we 
derived (using Eq.[TJ from the photometry. This difference is al- 
most certainly caused by effects of optical depth. 

8.2. Effects of optical depth r 

When analysing our data we tacitly assumed that the detected 
light originates, throughout the cycle, in the same gas and thus 
from the same atmospheric layer. Such assumptions have been 
made in almost all analyses of variable star data (also with 
the Baade-Wesselink method), but they need not be correct. 
Furthermore, Abt (1959) noted that the continuum light origi- 
nates in layers with t ^ 0.7 while the spectral lines rather are 
formed in layers with r ^ 0.3. 

During the cycle, compaction of the atmosphere and/or 
changes in level of ionization may lead to changes in t in that 
gas, thus to different layers being sampled in the photometry. 
This is also the case for spectral lines: changes in the density 
and temperature of a layer due to vertical motion will lead to 
changes in the local ionization balance and the excitation state of 
ions, thus to a varying strength of the relevant absorption lines. 
During the cycle a spectral line may therefore form in different 
gas layers. Furthermore, the strength of a line formed at greater 
depths is also influenced by the level of photon scattering filling 
in the spectral line. 

The phase lag between lines from metals and hydrogen dis- 
covered by van Hoof & Struve (1953) in a ySCepheid star has 
been attributed to similar optical depth effects. For RR Lyrae 
stars, this "van Hoof effect" was studied by, e.g., Mathias etal. 
(1995) and Chadid & Gillet (1998). Our sketchy spectral data 
(Fig. [8] Can and Hy) exhibit a phase lag too, which can easily 
be attributed to changes in the gas conditions that also cause the 
hysteresis in Teff and loggBj- 

In spectra of several RR Lyrae stars, line doubling has been 
reported near = (see, e.g., Sanford 1949). Considerable in- 
creases in line width may also occur near that phase (e.g., Feii 
lines; see, e.g., Chadid 2000). Spectral line doubling implies that 
there are gas layers in the same line of sight with the same ions 
but which are at different (radial) velocities. This doubling is 
seen mostly in the lower level Balmer series lines but also in 
other intrinsically strong lines such as Ca H&K. Thus, the spec- 
tral absorption in those lines comes from different geometric 
depths in the atmosphere. Note that when for the determination 



of radial velocities template spectra are used, any velocity differ- 
ences between ionic species are averaged out. 

8.3. Photometry and spectroscopy sample different layers 

Consider an outward moving dense layer of the stellar atmo- 
sphere. Upon its expansion and cooling, its continuum optical 
depth becomes smaller Thus, the level of r ^ 0.7 from which the 
continuum light emerges moves physically downward through 
the gas, leading to smaller photometrically derived R values. If 
this intrinsic r-level moves rapidly, the concomittant rapid de- 
crease in derived R is naturally interpreted as an extremely rapid 
downward velocity, even when the gas itself hardly has a vertical 
motion. The level from which the (metal) spectral lines emerge 
(near r ^ 0.3) will exhibit this effect later, i.e., only after the gas 
has been lifted further to larger radii (cooled and rarefied further 
to lower t). This happens later than the change in the gas levels 
releasing the continuum light. In this case, geometrically deeper 
levels may then be (spectroscopically) sampled suggesting a de- 
crease in R even if the gas itself still moves outward. We note one 
can see a difference in that sense between lines of Fe ii and Fe i 
(Chadid & Gillet 1998). Finally, the intrinsically strong lines, 
such as the lower Balmer series lines and a few other lines with 
large intrinsic absorption capability, will exhibit this effect yet 
later However, in deeper layers (at different radial velocity) the 
local Balmer absorption may already show up with its own ra- 
dial velocity (/l-shifted with respect to the Balmer line in higher 
layers), thus leading to line-doubling from lower levels of the 
Balmer series. 

The presented data do not allow us to asses the details of 
these optical depth effects. We emphasize that the parameters 
derived from the photometry always refer to the conditions in 
the layer with r ^ 0.7, the layer from which the measured con- 
tinuum light emanates. It may be that some of our interpretations 
need to be revised once more detailed measurements of the par- 
ticular gas layers would be available. 

The differences between our run of R (and of Vpui derived 
from it) and the Baade-Wesselink run of V^d (and of R derived 
from it) can be understood as being caused by the sampling of 
layers with different t. The Baade-Wesselink method samples 
layers with r smaller than r of the layers sampled with contin- 
uum photometry. It thus is unsurprising that the run of R as de- 
rived from these methods differs in important ways, all due to the 
physical condition (t) of the gas layer from which the measured 
radiation (be it spectral lines or the continuum) emerges. 

At the end of Sect. l7.2l it was noted that oscillations are visi- 
ble in Vpui with P/1. Perhaps these variations are also caused by 
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temporal small optical depth effects (Teff and loggBj) and do not 
represent fluctuations in gas velocity. 



9. The variations in log g(T,L, M) and log^Bj, and 
verifications of distance and mass 

The run of log g(T, L, M) was derived from Eq.[3] The actual sur- 
face gravities, log ^bj, were derived from b -y,ci. In Fig.fTOlwe 
plot three parameters: log g(T, L, M), log ^bj and the logarithmic 
ratio logigiT, L, M)/gsi). When log g{T, L, M) - log gsj < (or 
log{g(T, L, M)/gB]) < 1) the atmosphere is clearly compressed 
(collapsed). 

The data in Fig.[TO]indicate that in the quiet part of the cycle 
(0.4 < O < 0.8) the gravity logg(T,L,M) is, for several stars, 
not equal to log ^bj- However, in those quiet parts of the cycle, all 
parameters are apparently stable (there are most likely no optical 
depth effects at this point) and so we would expect log g(T, L, M) 
and log^Bj to be equal there. We can explore how one might 
adjust d and M to make log g{T, L, M) from Eq.[3]equal to log ^bj 
in that phase interval. For RR Lyrae stars, that phase (O ^ 0.6) 
has Teff ^ 6000 K, only a little higher than {T^s)q. 

The distances of our stars were defined as described in 
Sect. 121 They are taken from the literature and they are primarily 
determined from an adopted reference value of My- For some 
stars we calculated d from My (including effects of [Fe/H]). 

We recall from Sect. l4.5l (on the error budget) that the sum 
of the measurement errors in \og g{T, L, M) and loggBj may be 
as large as 0.08 dex. If these two gravities were to be equal 
and if we assumed the full margin of error of 0.08, then, e.g., 
\og g{T, L, M) for RR Gem should be reduced by ^0.22 dex, that 
of AS Cnc by ^0. 12, while that of TW Lyn should be larger by 
^0.07 dex. 

If the adopted My (or d) and M of our stars were wrong and 
had to be changed, this would produce (when using Eq.O the 
tabulated effects: 



assumed change 



effects 



Mv 


id) 


M 


\og g{T,L,M) 


L 


R 


-0.20 


+ 10% 


+ 10% 


-0.04 
+0.04 


+20% 


+ 10% 



We note again that changing the distance is for the stars observed 
only implicit; distances of RR Lyrae stars have (except for those 
in the papers referred to in Sect. l4.2l i not been determined by 
other means than through My- Thus only My or the mass M can 
be changed. 

For our stars a change in either My or M (or perhaps a com- 
bination thereof) is needed to make \og g{T,L,M) and loggBj 
match. The changes in My indicated below are within the range 
permissible by the non-uniqueness of My, as mentioned in 
Sect. 14. 51 For the stars shown in Fig.[TO]the changes are: 

- RR Gem: My 1 .0 mag brighter, M 50% smaller. The optimum 
would be to change My leading to L larger by a factor 2.5 (0.4 
dex in logL). Reducing the mass by 50% is no option, it then 
would be too low. A combination of changes might also fit. 

- TW Lyn: My 0.4 mag fainter, M 20% larger Changing My is 
not the right option because its L value would, in Fig.|6l be even 
further below the ZAHB. Hence M should be 20% higher at 0.85 
Mq. Note that (in Fig.|6]l the mean value of L can be made larger 
when adopting a yet larger mass. 

- AS Cnc: My 0.6 mag brighter, M 30% lower. Changes similar 
to but smaller than those for RR Gem would be needed. 

- SY Ari: no changes are needed within the errors. 

- SZ Gem: M smaUer by 0. 1 dex, to give M - 0.56 Mq. 
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Fig. 11. Three examples of Stromgren-photometry light curves 
of RR Lyrae stars near the possible "bump" in brightness. TW 
Lyn either does not show a bump or its bump is rather extended 
in phase; there is no or little effect on the colour indices. SY Ari 
does exhibit a bump in y but there are no effects in its colour 
indices. For AS Cnc, the colour indices u - b and u - v clearly 
change before y changes, indicating a change in the Balmer jump 
before a brightness change. The index b - y (representing tem- 
perature) changes somewhat later. 



- BH Aur: no changes are needed within the errors. 

- TZ Aur: to make the two gravities match, log g{T, L, M) should 
be ^0.3 dex lower. This would be achieved if M were consider- 
ably lower. Alternatively, My should (with M - 0.7 Mq) be so 
much brighter, that L is doubled. In both cases, TZ Aur would 
then be a well-evolved HB star 

We do not comment on the less well-observed stars (the bot- 
tom four in Table|2|. 



10. The lightcurve bump 

The stars in our programme were primarily selected from the 
class said to have lightcurve "bumps". These bumps in the run 
of brightness need not be very pronounced (for y-magnitudes 
see Fig.[T]l. One goal of the observations was to investigate 
whether these bumps in brightness affect the colours of the stars. 
Examples of the bump region for a few stars are given in Fig.fTTI 

First, we note that most of our stars with good light-curve 
coverage hardly show the light-curve bump in the Stromgren 
colour indices. Among these are TW Lyn, BH Aur, X CMi, TZ 
Aur, BR Tau, and RR Gem, which have either a prolonged low 
level bump or no clear colour-index bump at all. While SY Ari 
has a bump in y that has no effect in its colour indices, AS Cnc 
does show an effect. For the other stars, light-curve coverage was 
either poor or missing in this range of phases. 

In AS Cnc the brightening and dimming is clearly recog- 
nisable and takes place in a short timespan. We note that the 
colour indices u - b and u - v change before y brightens. Both 
indices become bluer meaning that the Balmer continuum radi- 
ation brightens without an increase in T^-s. In this part of the 
cycle, the atmosphere is cool and the Balmer continuum bright- 
ening could mean that the gas cools even further leading to an 
additional reduction in opacity based on the lower Balmer level 
excitation. Since T has not yet changed, the changes should be 
due to a reduction in gas density. Does extra light escape because 
of the lower r? At this point in the cycle, the atmosphere is quite 
extended and about to shrink. This happens at <1> ^ -0.1 (see 
Fig.[TT]and earlier ones). 

Bono & Stellingwerf (1994) attribute the bump feature to 
shockfronts. 
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1 1 . Further remarks 

11.1. Comparison with models 

Numerous models have been constructed for the pulsation and 
for various other detailed aspects of RR Lyrae (see, e.g., Fokin 
etal. 1999). A comparison with models is of limited value since 
most models are one-dimensional. More importantly, it is not 
always clear how theory has been transformed to observables. 
Do the theoretical predictions for V, R, spectral line strengths, 
and velocities all refer to information from the layer with the 
relevant r (thus to what one really would observe) or to a matter- 
defined gas layer? 

Parameters of RR Lyrae stars have been determined using 
dififerent kinds of data and with different modelling. We refer to 
the analysis of Sodor et al. (2009) for results based on the Baade- 
Wesselink method and its refinements as well as to additional 
literature. 

11.2. Blazhko effect 

A good portion of RR Lyrae stars shows cycle-to-cycle varia- 
tions in their light curves, the so-called Blazhko eff'ect (for ref- 
erences see, e.g., Jurcsik etal. 2009). Some of the stars of our 
sample are also known to be Blazhko stars but they exhibit these 
variations only at a moderate level. However, these variations 
may afi'ect the stellar parameters derived. The data we acquired 
cover to short a time span to address and assess possible cycle- 
to-cycle variations in the stellar parameters. 

1 1.3. Relation between stellar mass and kinematics 

Six of our stars were part of the study of kinematics of RR Lyrae 
stars (Maintz & de Boer 2005). Of these, five (CI And, AR Per, 
TZ Aur, RR Gem, and TW Lyn) are stars of the disk population 
(according to their kinematics), while SZ Gem is a star with halo 
kinematics. Halo RR Lyrae are understood to be older than disk 
RR Lyrae, so should (on average) have a lower mass. 

For TW Lyn, we found that its mass should be higher than the 
adopted reference value and instead be ^0.85 Mg, which chimes 
with it being a star of the disc. We note that Sodor et al. (2009) 
found a mass of ^0.82 Mq for RR Gem, our analysis indicated 
^0.7 Mq. For SZ Gem, we found a mass of ^0.56 M©, this low 
mass being in fine with the expected (larger age and) lower mass 
of halo stars. 

BH Aur and S Y Ari are found to correspond to the reference 
value of M = 0.7 M©, which is in line with their being disk stars. 

12. Conclusions 

Based on simultaneous Stromgren photometry and a compari- 
son of the photometric indices with a calibrated T^g, logg grid, 
accurate atmospheric parameters of RR Lyrae stars have been 
determined. Curves of phase related runs of log ggj (the gravity 
from the Balmer jump) and T^g were derived. One can then cal- 
culate the change in radius of the stars over the cycle. By includ- 
ing additionally obtained spectra, one arrives at the following 
description of the behaviour of RR Lyrae stars. 

The straightforward interpretation of the data derived for 
R is that at phase O ^ 0.9 the atmosphere begins to collapse 
soon reaching a high gas density (when the star has its greatest 
brightness). This manifests itself in a large Balmer jump (large 
log^Bj)- Within an interval A3) ^ 0.1 the atmosphere then ex- 
pands again. 



An alternative interpretation is that at phase 3> ^ 0.9 the op- 
tical depth of the extended outer layers has become so small that 
the level of r ^ 0.7 sampled in the photometry rushes through 
the gas inward mimicking a drastic reduction in radius. The con- 
traction taking place anyway then elevates the density in that 
level having r ^ 0.7, raising the optical depth there so that pho- 
tometry subsequently samples layers again higher up in the at- 
mosphere, suggesting a rapid expansion. The sudden increase in 
density manifests itself in lage values of loggBj- 

During all these changes the atmosphere exhibits oscillation 
ripples with a rhythm of P/7. 

A possible mismatch of the observed log^Bj and the cal- 
culated logg{T,L,M) in the quiet, descending part of the light 
curve can be used to assess the applicability of the mean param- 
eters adopted for the stars, i.e., the absolute brightness My (no 
geometric distances are known) and the mass M. One can then 
determine the individual values of My and M for a star 

Extending this kind of simultaneous Stromgren-photometry 
while avoiding multiplexing two stars would provide a much 
denser coverage of the light curves and thus a more accurate de- 
termination of the cycle variation in R. If spectra were then taken 
simultaneously, in which a range of spectral lines is included 
(lower and higher ionisation stages, strong and weak lines of the 
Balmer series) and of a nature to allow velocity determinations, 
one may hope to distinghuish the effects caused by the sampling 
of light from layers at difl'erent optical depth. 
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